Scheme to measure Majorana fermion lifetimes using a quantum dot 
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We propose a setup to measure the lifetime of the parity of a pair of Majorana bound states. The 
proposed experiment has one edge Majorana state tunnel coupled to a quantum dot, which in turn 
is coupled to a metallic electrode. When the Majorana Fermions overlap, even a small relaxation 
rate qualitatively changes the non-linear transport spectrum, and for strong overlap the lifetime can 
be read off directly from the height of a current peak. This is important for the usage of Majorana 
Fermions as a platform for topological quantum computing, where the parity relaxation is a limiting 
factor. 

PACS numbers: 74.25.F-, 85.35. Gv, 74.45. +c, 74.78.Na, 



Topological superconductors are currently attracting 
massive interest, partly due to the Majorana bound 
states (MBS) which form at edges or vortices of such sys- 
tems. Majorana Fermions are non-Abelian anyons 
meaning that particle exchanges are non-trivial opera- 
tions which in general do not commute and can be used to 
perform quantum computational operations, called topo- 
logical quantum computing 0]. Non-Abelian statistics 
has been predicted for the quasi-particle excitations of 
the v = 5/2 fractional quantum Hall state J3|, and ap- 
pear as a consequence of p-wave type pairing jj] , and such 
superconductors should host MBS in vortices. Recently, 
it was realized that p-wave like pairing may also occur 
in topological insulators and even in ordinary semi- 
conductors with strong spin-orbit coupling w hen 
brought into proximity with an s-wave superconductor. 

A particularly simple system is a one-dimensional 
semiconducting wire with strong spin-orbit coupling, 
brought into proximity with an s-wave superconduc- 
tor Ql, Hp. The original idea of MBS in wires is much 
older the more recent proposals being possible ex- 
perimental realizations of that model system. Particle 
exchanges (braiding) could be accomplished by crossing 
two or more wires [12l Il3j. However, the first step would 
be to verify the existence of MBS, which could e.g., be 
achieved by tunnel spectroscopy 0, [t| . The presence of 
a MBS gives rise to a characteristic zero-bias conductance 
peak, while a more complicated peak-structure arises if 
many MBS are coupled to each other (l6| . 

The main advantage of topological quantum comput- 
ing is the robustness against decoherence: the compu- 
tational basis consists of pairs of MBS which are spa- 
tially separated and the state normally cannot be mod- 
ified by perturbations which do not couple simultane- 
ously to more than one MBS 0. However, perturba- 
tions which change the parity degree of freedom of the 
superconductor can change the state of the Majorana 
system and thus lead to decoherence Q . The presence of 
such p arity -changing processes, called quasiparticle poi- 
soning |l7l [l8j , is in fact a well-known problem in super- 
conducting charge qubits. In order to determine whether 
Majorana bound states in different topological materials 
are actually suitable for topological quantum computing, 
it is therefore crucial to be able to measure the lifetime 




FIG. 1: (Color online) Sketch of the proposed setup. An s- 
wave superconductor induces superconductivity in a semicon- 
ducting wire with strong spin-orbit coupling. A magnetic field 
can then induce a topological phase where MBS are formed 
at the ends of the wire. A quantum dot is formed e.g., by 
separating a short segment from the main wire by a tunnel 
barrier, and it is tunnel coupled to a metallic normal lead. 



of the parity degree of freedom. Here we show how this 
can be done in an experimentally rather simple transport 
setup, which does not require braiding or interferometry. 

The parity lifetime cannot be measured in a normal 
metal-MBS tunnel junction as studied in Refs. [T34l6| . 
and we consider instead a setup including a quantum 
dot between the normal metal and the topological su- 
perconductor, sec Fig. [TJ We focus here on a wire-type 
topological superconductor, which seems experimentally 
most attractive, but the main conclusions are indepen- 
dent of such details. Systems of several coupled quan- 
tum dots and MBS have been suggested as probes of 
the non-locality of Majorana Fermions and very re- 
cently as an alternative way to perform non-Abelian rota- 
tions within the degenerate Majorana groundstate man- 
ifold (20| . However, the non-equilibrium transport prop- 
erties of a normal metal-quantum dot-MBS junction 
have to our knowledge not been investigated previously 
and therefore the possibility to probe the parity lifetime 
in such a setup has also not been discovered. As we 
will show, even a small coupling between the two MBS 
strongly suppresses the current. A finite parity relax- 
ation rate partially restores the current but qualitatively 
changes the transport spectrum, allowing a finite lifetime 
to be detected, and in addition its precise value can be 
measured if the MBS coupling or tunnel rates can be 
controlled. 

We solve exactly the problem of a strongly interacting 
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quantum dot coupled to a MBS and treat the coupling to 
the normal metal perturbatively, yielding a set of master 
equations for the state of the combined dot-MBS system. 
Prom these equations we calculate the non-linear current 
as a function of the bias voltage applied to the normal 
lead and of the gate voltage which controls the quantum 
dot energy. In addition, in the supporting information 
we provide an exact solution of the transport problem 
for the case of a non-interacting quantum dot [21] . 

The basic properties of MBS in semiconducting sys- 
tems with induced superconductivity have been de- 
scribed in e.g., Refs. @, 0,111 US- For suitable parameters, 
MBS are formed at the ends of the wire, see Fig.[TJ These 
are zero-energy solutions to the Boguliubov de Genne 
equations and are described by the operators 

7 m = / <te$^(Som(«)*a(«) + g* m (z)*t(z)) , (1) 

where ^\(z) creates an electron with spin projection 
(7 4- at position z in the wire. The exact form of the 
envelope functions j, m (z), describing the spatial form of 
the MBS, depends on the details of the wire, but g a \ and 
g a 2 have their main weight close to opposite ends and 
decay exponentially inside the wire. Clearly, the MBS 
operators satisfy 7 m = 7 m and we assume them to be 
normalized, such that 7172 = —7271, 7m = 1- One end 
of the wire is tunnel coupled to a quantum dot. The 
coupled dot-MBS system is described by 

H = H D + ^7i72 + £ {Kd a - A* 4) 71, (2) 

where Ht, = e <J n <j + Un^n^ describes the dot, n a = 
d^da is the number operator, and U is the Coulomb 
charging energy for electrons on the dot. The low-energy 
Hamiltonian describing the two MBS was discussed in 
Ref. [ill and has been used in a transport setup in e.g., 
Refs. [Il 03-G1 El; the coupling of the two MBS is 
given by £, and A CT is the coupling between electrons 
on the dot and the MBS at the corresponding end of 
the wire. To be specific, we assume the wire to have 
a Rashba-type spin-orbit coupling along the x-direction, 
and apply a magnetic field B along z, leading to a Zee- 
man splitting on the dot, = £f — ej,. For definiteness, 
we take the MBS wave function at the dot side of the 
wire to be oc i z ) + ^i y ( z ) ( see e -S-, 01 )> where i y 
means spin along the negative y-axis, and therefore cou- 
ples only to the y-component of the dot spin, equal in 
amplitude but different in phase for the two dot spin di- 
rections, A-f = A, Xi = —iX. The metallic normal lead is 
described by Hjf = J2ka efcC lo- Cfccr an< ^ ^ C0U Pl es to the 
dot via Ht = Y^ka^^^ka + h.c, where c' k creates an 
electron in the normal lead with energy e^, and where 
is the amplitude for dotonormal lead tunneling. 

It is useful to switch to a representation where the two 
Majorana Fcrmions are combined to form one ordinary 
Fermion: 71 = / + /t , 72 = — /), where p creates 



a Fermion and p f = 0, 1 counts the occupation of the 
corresponding state. The Hamiltonian ^ now becomes 

Ha = H D +^pf-^j 

+ (Aada/ f + Kdaf + h.c) , (3) 

When A = 0, the eigenstates of ^ are given by \nd,nf), 
where = 0,^,1,2 and nj = 0,1 describe the states 
of the dot and the MBS system, respectively. Also when 
A / 0, the Hamiltonian is block-diagonal in this basis, 

with an even block H^\ acting on |00), \<jl), [20) , and 

an odd block Hq , acting on |ct0), |01), |21). Fermion 
number is not conserved by ([3]), but, because of the block- 
structure, the parity of the total Fermion number of the 

(e) 

dot-MBS system is conserved. When £ = 0, Hq ' and 

H^ are identical, £ ^ breaks this parity symmetry. 
Diagonalizing ^ yields odd/even parity eigenstates 

\oi/e l ) = T, nd n f eo/e a nIn f \n d n f ). The tunnel Hamil- 
tonian, Ht, changes the dot electron number by ±1 
and thus connects the even and odd parity sections: 
(ei\H T \oj) = J2ka K ka' Aa + h-c, where the many-body 
tunnel matrix elements are given by 

K tT = ** £ a nan f a°n' d n' f { n d^f\da\n d n f ). (4) 
ndnfQe 
n' d n' f £o 

In the presence of strong electron-electron interactions, 
the tunneling between the normal lead and the dot-MBS 
system cannot be solved exactly (in the limit [7 = 
we provide an exact solution in the supporting infor- 
mation [2l|). To leading order in Ht, the problem can 
be formulated in terms of master equations [22l |, which 
we modify to account for the lack of particle conserva- 
tion. The occupations P a of the dot-MBS eigenstates 
| a) = |e,-), \oi) are calculated from 

= ( W ™' P »' - W a , a P a ) , £ P a = 1, (5) 

a' a 

W aa> = J2{ T f F ( E aa')+r a J a [l~F(E a , a )}},(6) 
a 

where E aa i = E a — E a i , with E a the energy of eigenstate 
\a), and F(E) = \/[e^ E -^)l T + 1) is the Fermi func- 
tion of the normal lead with electron temperature T and 
chemical potential /ijv = V. (We set kg = e = h = 1.) 
The tunnel couplings arc T® a = 2irpN\K% a | 2 , with 
being the density of states in the normal lead. Note 
that due to the lack of particle conservation, any pair 
of states a, a' are connected both by processes removing 
and adding an electron to the normal lead, described by 
the first and second term in Eq. ([6]) respectively. We here 
assumed energy independent density of states and tun- 
nel amplitude, but the eigenstates a, a' depend on the dot 
level position and the effective tunnel couplings T^ a vary 
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between and 2ivpN\t\ 2 — T. The physics described by 
Eq. (J5J) is easily understood. The occupation probabil- 
ity of state a is given by the sum of all tunnel processes 
starting from any state a' and ending with occupation 
of state a (first term) , minus all processes depopulating 
state a (second term). The occupations are normalized 
to one. Note that the rates in Eq. describe electron 
tunneling only, rates related to parity relaxation can be 
added to W aa ' as described below. The particle current 
flowing into the dot from the normal lead is given by 
^ = Saa' Woa'-Pa'i where the current rate matrix W* a , 
is similar to (O, but with a minus sign on the second 
term, corresponding to electrons tunneling into the nor- 
mal lead. 

The restriction to lowest order perturbation theory is 
valid when T»T. In addition, we neglected off-diagonal 
elements of the density matrix, valid when all eigenstates 
are well separated on a scale set by the tunnel broaden- 
ing, satisfied when A ^> T. Neglecting all states in the su- 
perconductor except the MBS is valid when the induced 
superconducting gap A is larger than the other energy 
scales. Moreover, U can be several meV in small quan- 
tum dots formed in semiconducting wires [H, [24| , which 
is likely > A, and the same holds for the level-spacing, 
motivating the single-orbital model used here. 

Next, we present solutions of the master equations for 
specific choices of parameters, and with the above en- 
ergy scales and limitations in mind, we take U = oo 
and T = A = 10r = An/100, with T < 100 mK. The 
voltage-dependence of the energy cost for adding an elec- 
tron to the dot can be included in the orbital energies, 
(e± + 6f)/2 = eo — &Vg + [3V, where a (gate coupling) 
and (3 depend on the capacitances associated with the 
dot-superconductor and dot-normal metal tunnel junc- 
tions. We take (3 = 1/2 and eo = (other choices affect 
only the slope and absolute position of resonances, re- 
spectively). The differential conductance, dl/dV, plot- 
ted on color scale as function of V and V g ("stability 
diagram"), is shown in Fig. [5] for different values of A 
and £. 

To develop an understanding of the physics underlay- 
ing the stability diagrams, we consider first the case with- 
out parity relaxation, described by Eqs. flS) and (JS). We 
start with the simplest case of £ = 0, sec Fig.[2fa). A sig- 
nificant current can only flow when one of the dot spin 
states is energetically aligned with the MBS, allowing 
electrons to tunnel resonantly between the dot and su- 
perconductor, resulting in a current peak of width ~ A 
and height / max =T/2 and therefore in narrowly spaced 
positive and negative conductance lines. Thus, the sharp 
nature of the MBS gives rise to negative differential re- 
sistance in transport through the dot. The peak asso- 
ciated with a =t terminates at large positive V, when 
an electron can tunnel from the normal lead into the 
(7 =1 state, blocking transport through the a =t chan- 
nel (via the Coulomb charging energy) since it cannot 
tunnel out again into the normal lead (due to energy 
conservation) or into the MBS (due to the energy mis- 




(e) |oo) t* [ffQ) |^o> x > |oi) 




FIG. 2: (Color online) Stability diagrams: dl /dV normalized 
by 0.02e 2 //i, plotted on color scale as function of V and V g , 
with (a): £ = 0, (b): £ = 10T, (c): £ = 10T and finite parity 
relaxation, A = 0.04I\ Note the different scales in (b) and 
(c) showing the dramatic effect of the finite parity lifetime 
and the characteristic shift in gate voltage position between 
the positive and negative bias peaks in (c). (d): I(V g ) at 
V = 75T with parameters as in (c), but increasing f/A = 
0, 2.5, 5, 10, 20, 30 (increasing as indicated by the arrow). For 
can be read off directly from the height of 
the shifted current peaks, (e): Sketch of resonant transport 
processes for e = £ when V > in (c). Similar resonant 
processes exist for e = — £ when V < 0. 



match). Similarly, the a =4- peak vanishes above a neg- 
ative threshold V, when a =\ falls below /xjv- Thus, 
properly accounting for Coulomb blockade is essential 
to understand the qualitative conductance features. In 
Fig. |3£b), a finite £ introduces an energy splitting be- 
tween the even and odd parity sectors, suppressing the 
current for \V\ < £. Even for \V\ > £ the conductance 
is significantly lower compared to £ = [note the dif- 
ferent scales in Figs. (2ja) and (b)]. To understand why, 
consider e.g., electrons being transported from the nor- 
mal lead to the superconductor (V > 0) by sequential 
tunneling through the dot, which in the basis of the un- 
perturbed number states \ndUf) involves the processes: 
|00) "7* |er0) ~? |01) |oT) |00). Because we con- 
sider here the case without parity relaxation in the super- 
conductor, any process changing the parity of the dot- 
MBS system must involve tunneling to/ from the normal 
lead, and therefore two electrons are transferred before 
the system returns to its initial state. When the even- 
odd degeneracy is split by £ ^ 0, both dot-MBS tunnel 
processes cannot be resonant at the same voltage, sup- 
pressing the conductance by a factor oc A 2 /£ 2 . Instead, 
horizontal conductance lines appear, corresponding to in- 
elastic cotunneling [2|| through the dot, exciting the par- 
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ity degree of freedom, i.e., a direct processes |00) ~p\ 1 1 ) , 
energetically allowed for \V\ > £ and involving only vir- 
tual occupation of the intermediate state |<r0). Note that 
all tunnel process oc TX n are included in our master equa- 
tions, since the dot-MBS coupling is treated exactly. In- 
terestingly, due to the sharp MBS, inelastic cotunneling 
gives rise to conductance peaks, rather than steps. Simi- 
larly, horizontal lines at V = ± correspond to inelastic 
cotunneling exciting the dot, e.g., | t 0) at | i 0). 

Next, we introduce a finite relaxation of the Fermion 
parity in the superconductor, caused by quasiparticle poi- 
soning [It], [H| j an d show that this has a striking effect on 
the stability diagrams. Rather than starting from a mi- 
croscopic model of the quasiparticle generation, we focus 
on the generic effects and assume that transitions be- 
tween dot-MBS cigenstatcs \a'} —> \a) are induced with 
rate A aa < if E a , > E a and rate A aa >exp[(E a * - E a )/T] 
otherwise, where A aa i is obtained by assuming a rate 
A for transitions \ndUf) — > \ridfif) (1 = 0,0 = 1), 
and then transforming to the dot-MBS eigenbasis, cf., 
Eq. ([!]). These rates are then added to the rate matrix 
in Eq. ©. As is seen by comparing Fig. [^c) and (b), 
which differ only by a small such relaxation rate, this 
has a dramatic effect on the conductance (note the dif- 
ferent scales). To understand the increased conductance, 
consider transport at V > 0, which for A ^ is dom- 
inated by the processes: 1 00) "J* |er0) 1 1 ) ~^ |00), 
see sketch in Fig. [He). Hence, parity relaxation allows 
the total dot-MBS system to return to its initial state 
after transferring only a single electron, and this elec- 
tron tunnels through the system fully resonantly when 
e a = £. Similarly, at V < transport is dominated by: 
1 00} ~x | erf) \a0) ~p* 1 00) , which is fully resonant 
at e a = — £. Thus, parity relaxation leads to resonant 
transport even for finite £, but the resonance condition is 
different for positive and negative bias and therefore the 
peaks are shifted by 2£. 

The shifted peak behavior is observed when the parity 
relaxation rate exceeds the suppressed transport rates 



discussed in connection with Fig. [2][b), i.e., when A > 
T\ 2 /£ 2 . If, in addition, r> A, parity relaxation is the 
slowe st p rocess in Fig. [^e) and the peak current is equal 
to A [26| . Experimentally, £ could be changed by moving 
the MBS corresponding to 72 with "keyboard" gates as 
suggested in Ref. [l3|]. I(V g ) curves for increasing £ are 
shown in Fig. [21(d), allowing the parity relaxation to be 
directly read off from the height of the current peaks at 
large £. Alternatively, one could change instead T and/or 
A via additional gates. Note that A can be directly mea- 
sured without knowing the precise values of e.g., T and 
A. We emphasize that even at finite £, the two-Majorana 
state described by 71,2 or fi,f can only relax by also 
changing the parity of the superconductor. Therefore, 
even though our proposed measurement has to be done 
at finite £, the relaxation rate measured is the parity re- 
laxation rate, which is arguably independent of £. 

In conclusion, we have suggested a relatively simple 
setup in which to measure the lifetime of the parity of 
a Majorana system, being of crucial importance in topo- 
logical quantum computing schemes 0. Our suggestion 
does not involve (experimentally problematic) interfcr- 
ometry or braiding operations. Thus, one can test differ- 
ent realizations of topological superconductors and dif- 
ferent material choices, selecting the most suitable ones, 
before going through the difficult process of actually con- 
structing a topological quantum computer. The setup 
we propose is a normal metal-quantum dot-topological 
superconductor junction. We have calculated the non- 
linear conductance of this system as a function of the 
applied gate and bias voltages within a master equation 
approach and shown that the lifetime of the parity of the 
Majorana system can be directly read off from the height 
of the current peaks at finite Majorana coupling. 

We thank Harvard University, where part of this work 
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by The Danish Council for Independent Research | Nat- 
ural Sciences and by Microsoft Corporation Project Q. 
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where Fq(lj) is the Fermi function at V = 0, G dd 
is the Fourier transform of the retarded dot Green's 
function for a =t electrons, and F dd (t — t') = 
—i0(t — t')({df(i),d^(t')}) is the corresponding anoma- 
lous Green's function. (We here only consider A = 0.) 
Next, we take the time derivative of G dd (t — if) [and of 
all Green's functions appearing when doing so, including 
F^it — t')], Fourier transform, and solve the resulting 
equations of motion, which gives 



w-e t + fr-2|A| 2 M(o;) 1 + 2|A| 2 M (w) 



(8) 
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FIG. 3: (Color online) (a) dl/dV normalized by 2e 2 /h, plot- 
ted on color scale as function of V and V g , with £ = T = 2A, 
T = 0. We take the voltage coupling j3 = 1/2 as in Fig. 2 in 
the main paper, but choose the gate voltage scale such that 
6f = when V g = V = 0. (b) dl/dV as function of V at 
V g — 0. The parameters are the same as in (a), except for £ 
which is varied: £ = (red solid line), £ — T/2 (green dashed 
line), £ = F (blue dotted line), £ = 3F/2 (purple fine dotted 
line) . 

We here consider the Zecman splitting to be the largest 
energy scale, As 3> \V\,T, T,A, in which case transport 
through the ground state of the dot can be described by a 
non-interacting model and solved exactly using a Green's 
function approach. The current can be written as 

I = -r/^ [2ImG« +T(\G^-\F^)] 
x [F(ui)-F (u})], (7) 



-2{\*) 2 M{u) 

(w + e t + |r) (w - e t + |r) - 4|A| 2 (w + |r) M (w) ' 

(9) 



where M(u) = l/(w - £ 2 /l<j), and M(w) = M(lu)/[lu + 
e t + §r-2|A| 2 M(w)]. 

When A, T ^> T, the results agree fully with the master 
equation results for B = oo and we now focus instead on 
the limit T = 0, A < V. Figure El^a) shows a stability 
diagram similar to those in Fig. 2 in the main paper, but 
for a smaller voltage window close to the e-j- = crossing. 

Figure |31[b) shows dl/dV as function of V at V g = 0, 
for varying values of £. The zero-bias conductance, G(0), 
can be found analytically from (0. When £ = 0, we al- 
ways have G(0) = 2e 2 /h, independent of V g , although 
the width of the peak decreases when the dot level is 
brought out of resonance. Conversely, for any £ ^ 0, we 
find G(0) = 0, independent of V g . This is similar to the 
linear conductance of a normal metal-MBS tunnel junc- 
tion, sec Ref. [16] of the main paper. At finite £, peaks 
instead appear at finite voltage V ~ ±£, as discussed in 
connection with the master equation results in the main 
paper, and close to V g — an interesting double peak 
structure is found. 



